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Abstract 

State-of-the-art ab initio techniques have been applied to compute the potential energy curves 
of the (BaRb) + molecular ion in the Born-Oppenheimer approximation for the singlet and triplet 
states dissociating into the ground state 1 S Rb + ion and the Ba atom in the ground *S state, the 
lowest singlet or triplet D excited states, and for the singlet and triplet states dissociating into the 
ground state 2 S Rb atom and the ground state 2 S Ba + ion. The ground state potential energy 
was obtained with the coupled cluster method restricted to single, double, and nonperturbative 
triple excitations. The first triplet states in the S, II, and A symmetries were computed with 
the restricted open-shell coupled cluster method restricted to single, double, and nonperturbative 
triple excitations. All other excited state potential energy curves were computed using the equa- 
tion of motion approach within the coupled-cluster singles, doubles, and linear triples framework. 
The long-range coefficients describing the electrostatic, induction, and dispersion interactions at 
large interatomic distances are also reported. The electric transition dipole moments governing the 
X X S — > 1 S, 1 I1 have been obtained as the first residue of the polarization propagator computed 
with the linear response coupled-cluster method restricted to single and double excitations. Nona- 
diabatic radial and angular coupling matrix elements, as well as the spin-orbit coupling matrix 
elements have been evaluated using the multireference configuration interaction method restricted 
to single and double excitations with a large active space. With these couplings, the spin-orbit 
coupled (relativistic) potential energy curves for the + and 1 states relevant for the running ex- 
periments have been obtained. Finally, relativistic transition moments and nonadiabatic coupling 
matrix elements were obtained from the nonrelativistic results and spin-orbit eigenvectors. The 
electronic structure input has been employed in the single channel scattering calculations of the 
collisional cross sections between the Ba + ion and Rb atom. Both nonrelativistic and relativistic 
potentials were used in these calculations. Our results show that the inelastic cross section corre- 
sponding to the charge transfer from the Rb atom to the Ba + ion is much smaller than the elastic 
one over a wide range of energies up to 1 mK. This suggests that sympathetic cooling of the Ba + 
ion by collisions with ultracold Rb atoms should be possible. 
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I. INTRODUCTION 



Nowadays an increasing number of experimental groups worldwide start to work with 
hybrid systems involving cold or ultracold trapped atoms and ions jl-5]. Apart from the 
fundamental interest in the physics of atom-ion collisions in the quantum regime 6|-[8| , these 
ne w sy stems are very attractive from the point of view of quantum information processing 
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101 ] . studying many-body effects of ion impurities [llj, or creation of the molecular ions 



12|. One type of ongoing experiments is focused on studying the co.lis.ona. processes in 
cold clouds of atoms and ions, like Yb with Yb + and Na with Ca + [l|, [2[ stored in dual 
hybrid charged-neutral traps at mK temperatures. Other experiments, in contrast, study 
the dynamics of single ions like Ba + or Yb + trapped in RF potentials and immersed in Bose- 
Einstein condensates In such experiments the ultracold cloud of atoms is prepared 

using standard cooling and trapping techniques, while the ion is laser cooled separately in 
RF trap and later overlapped with the Bose-condensed atoms. The atom-ion collisions can 
lead to further sympathetic cooling of the ion, but the net cooling effect depends on the 
interplay of (i) the two-body collisional properties, (ii) the micromotion of the ion due to 
the time-dependent RF potential, (Hi) the collective phenomena resulting from the coherent 
properties of the condensate additionally modified by the presence of an ionic impurity. In 
this paper we perform a first step towards understanding the sympathetic cooling process of 
Ba + ion with Rb atoms, by calculating highly accurate molecular potentials and determing 
the single channel elastic, spin-exchange and charge-transfer collision rates. 

In a recent paper Makarov et al. 13) performed ab initio and dynamical calculations on 
the (CaNa) + molecular ion. The results of these calculations suggest that the milliKelvin 
regime of collisional cooling of calcium ions by sodium atoms is very favorable, with the rate 
coefficient for charge transfer from the Na atom to the Ca + ion several orders of magnitude 
smaller than the rate for elastic and spin-exchange collisions. This strongly suggests that 
sympathetic cooling of ions by collisions with ultracold atoms should be possible. This 
system was further studied in the ultracold regime in Ref. [8j] within the multichannel 
quantum defect formalism. In the present paper we investigate the possibility of sympathetic 
cooling for yet another system of experimental interest [5]: Ba + ions cooled by collisions with 
ultracold Rb atoms. 



Theoretical modeling of collisions in the ultracold regime requires a lot of care 14J4l8|. 
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First of all, the electronic states involved in the dynamics must be computed with the 
state-of-the-art methods of quantum chemistry. In particular, these methods should be 
size-consistent in order to ensure a proper dissociation of the molecular state, and must 
sufficiently account for the electronic correlation. Moreover, any information on the long- 
range asymptotics of the potentials is very important. Finally, all couplings between the 
molecular states, both those resulting from the spin-orbit interaction and from the nonadi- 
abatic effects, should be considered. Having the electronic structure results at hand, exact 
quantum-dynamical calculations should be performed to get the cross sections and collisional 
rates. 



In a recent series of papers Knecht et al. 19|, |20) reported nonrelativistic, scalar relativis 



tic, and fully relativistic ab initio potential energy curves for the (BaRb) + molecular ion. 
Unfortunately, the approach adopted in these papers is not size-consistent, so the results 
at large internuclear distances may not be accurate enough. Moreover, excitations in the 
wave function beyond the doubles levels were included only for the simplest single reference 
closed-shell and high-spin states. Also, the transition moments between states, necessary 
to model the radiative charge transfer process from the Rb atom to the Ba + ion, were not 
computed. Finally, the nonadiabatic coupling matrix elements between different molecular 
states were not considered. We believe that such a level of ab initio calculations is not suffi- 
cient to properly model ultracold collisions in the (BaRb) + molecular ion. Therefore, in the 
present paper we report state-of-the-art ab initio potential energy curves for the (BaRb) + 
molecular ion in the Born-Oppenheimer approximation for the singlet and triplet states dis- 
sociating into the ground state X S Rb + ion and the Ba atom in the ground X S state, the 
lowest singlet or triplet D excited states, and for the singlet and triplet states dissociating 
into the ground state 2 S Rb atom and the ground state 2 S Ba + ion, electric transition dipole 
moments, nonadiabatic and spin-orbit coupling matrix elements. Except for the spin-orbit 
coupling and some nonadiabatic coupling matrix elements, all the results were obtained with 
size-consitent methods based on the coupled cluster ansatz including triple excitations. The 
nonrelativistic results are transformed to the relativistic basis, thus allowing to judge the 
importance of the relativistic effects, and to use our data in quantum dynamical calculations 
both within the Hund's case (a) and (c). Finally, we report single channel calculations of the 
elastic, spin-flip, and charge transfer cross sections. This will allow us to obtain the ratio of 
inelastic to elastic cross sections at various temperatures, and thus give a first estimate of 
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the efficiency of the sympathetic cooling of the barium ion by collisions ultracold rubidium 
atoms. 

The plan of this paper is as follows. In Sec. [TT] we introduce the theoretical models used 
in our calculations. We start this section with a description of the methods used in ab initio 
calculations of the Born-Oppenheimer potential energy curves and electric transition dipole 
moments for (BaRb) + . We continue with the calculations of the nonadiabatic and spin-orbit 
coupling matrix elements, and of the relativistic potentials. The choice of fixing the long- 
range coefficients at their ab initio values is also addressed. The remaining part of this section 
is devoted to the second step of the Born-Oppenheimer approximation, i.e. to dynamical 
calculations of the elastic, spin-flip, and charge transfer scattering cross sections. Numerical 
results are presented and discussed in Sec. IHIl We start this section with the discussion of 
the ground and excited states potentials. Next we turn to the nonadiabatic coupling matrix 
elements and electric transition dipole moments. The effects of the spin-orbit coupling on 
the potentials, transition moments, and nonadiabatic coupling matrix elements are also 



addressed. W 
initio data 19 



renever possible, our theoretical results are compared with the available ab 



20J]. Once the results of the electronic structure calculations are presented 
and discussed we turn to the problem of producing cold barium ions. We present the results 
for the elastic, spin-flip, and charge transfer cross sections for both the nonrelativistic and 
relativistic potentials, and discuss the efficiency of the sympathetic cooling leading to cold 
barium ions. Finally, in Sec. IIVI we conclude our paper. 

II. COMPUTATIONAL DETAILS 
A. Electronic structure calculations 

When dealing with collisions at ultra-low temperatures the accuracy of the potential in 
the long range is very important. Therefore, the methods used in the calculations should be 
size-consistent in order to ensure a proper dissociation of the electronic state, and a proper 
long-range asymptotics of the potential. In the present paper we adopt the computational 
scheme successfully applied to the ground and excited states of the calcium dimer 



14 



181. 



The potential energy curves for the ground and excited states of the (BaRb) + molecular ion 
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have been obtained by a supermolecule method: 

V 2S+1 ^(R)=E™-El M -E* M (1) 

where -E^b denotes the energy of the dimer computed using the supermolecule method SM, 
and £f M , X=A or B, is the energy of the atom X. For the ground state potential we used 
the coupled cluster method restricted to single, double, and noniterative triple excitations, 
CCSD(T). For the first triplet states of the E, II, and A symmetries we employed the 
restricted open-shell coupled cluster method restricted to single, double, and noniterative 
triple excitations, RCCSD(T). Calculations on all other excited states employed the linear 
response theory (equation of motion) within the coupled-cluster singles, doubles, and linear 



triples (LRCC3) framework [2l|, |22j. Note that the second and higher triplet states and all 
excited singlet states are open-shell systems that cannot be described with a single high-spin 
reference function, so one has to resort to methods especially designed to describe open-shell 



situations 
program 



231 ] . The CCSD(T) and LRCC3 calculations were performed with the dalton 



24j . while RCCSD(T) calculations were done with the molpro suite of codes (25]. 
It is interesting to note at this point that even though (BaRb) + is effectively a two-electron 
system, triple excitations are very important. A more detailed discussion of this point will 
be presented in sec. IHIl In principle we could use the LRCC3 method to obtain all the 
excited state potentials in any symmetry. However, the LRCC3 method is computationally 
more expensive than the RCCSD(T) approach, so we decided to use the latter when possible. 
However, we have checked for a few points that the RCCSD(T) and LRCC3 results for the 
(1) 3 S, (l) 3 n, and (1) 3 A are very close. 

The long-range asymptotics of the potentials is of primary importance for cold collisions. 
Therefore, for each state we have computed the leading long-range coefficients describing the 
electrostatic, induction, and dispersion interactions. For the ground X X S state dissociating 
into Rb + ( 1 S)+Ba( 1 S) and for the first excited singlet and triplet states dissociating into the 
Ba + ( 2 S)+Rb( 2 S) the leading long-range asymptotics at large internuclear distances R reads: 

si'md siind /^ydisp 
2S+i|A| , . Q 4 Q 6 Og / 9 n 

where the lone-range coefficients Q nd , Q nd , and C 6 disp are given by the standard expres- 



sions (see, e.g. Refs. [28|, |29j) that can be derived from the multipole expansion of the 



interatomic interaction operator. The long-range induction and dispersion coefficients were 



computed with the recently introduced explicitly connected representation o 



;he expecta- 



tion value and polarization propagator within the coupled cluster method [30|, |3l|, and the 
best approximation XCCSD4 proposed by Korona and collaborators 32)]. For the singlet 
and triplet states dissociating into Ba + ( 2 S) ion and Rb( 2 S) ion the induction coefficients 
were obtained from finite- field RCCSD(T) calculations, while the dispersion coefficient from 
the sum-over-state expression with the transition moments and excitation energies computed 
with the multireference configuration interaction method limited to single and double exci- 
tations (MRCI). Specifically, the transition moments and excitation energies of the Ba + ( 2 S) 
ion were obtained in this way, while the Rb polarizability at imaginary frequency was taken 
from highly accurate relativistic calculations from the group of Derevianko [3^] . 

For the molecular states of the (BaRb) + ion dissociating into Ba( 1 / 3 D)+Rb + ( 1 S) the 
long-range asymptotics of the potentials is slightly more complicated and reads: 

/^<elst /^i'md /^<elst z^i'md /^rdisp 
I/ 2S+1 |A|/n\ °4 1 °5 °6 ^6 , ( o\ 

where the new terms appearing in the expression above describe the long-range charge- 
quadrupole (C| lst ) and charge-hexadecapole (C| lst ) interactions. The mathematical expres- 
sions for the coefficients of Eq. ([3]) are given by: 



Cf st = (-1) 2+A ( 1 / 3 D||g 2 || 1 / 3 D), (4) 




2 V 6 

Ct = (-ir A ( 2 4 H ( 1/3 D|iQ 4 |i 1/3 D), (6) 
\-A A J 

Cind ^ i i 

6 — 2 zz ' ZZl ^ ' 

Cf* = ^l°° «o Rb+ (^) (teoiiv) + 3^y^a 2 M) du. (8) 

In these equations the expression in round brackets is a 3j symbol, ( 1//3 D| \Qi\ | 1//3 D), (7 = 2 
or 4) denotes the reduced matrix element of the quadrupole and hexadecapole moment, 
respectively, ao and are the scalar and tensor components of the electric dipole polariz- 
ability tensor of the Ba atom in the 1//3 D state 33J], while C ZZ)ZZ is the z component of the 



quadrupole polarizability of Ba( 1 / 3 D). Finally, a^ b+ is the polarizability of the rubidium ion 
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Rb + ( 1 S). The electrostatic coefficients C| lst and Cf st for the E, II, and A states are not 
independent, and are connected one to the other by the following relations: 

C 3 olst (E) = -Cf st {A) = 2Cf rt (II), Cf s = 6Cf *(A) = -^Cf st {U). (9) 

Note parenthetically that all coefficients that lead to attractive interactions (induction and 
dispersion terms) are assumed to be positive, while the electrostatic constants may result 
in both attractive and repulsive interactions, so they enter with their true sign. The values 
of the quadrupole moments, and scalar and tensor dipole polarizabilities, as well as the 
components of the quadrupole polarizability were obtained from finite-field LRCC3 calcu- 
lations on the 1/,3 D state of the atom. The hexadecapole moment of the singlet state was 
obtained as an excited state expectation value within the linear response CCSD formalism of 
Christiansen et al. 26], while for the triplet state from the MRCI calculations. The vector 



component of the dipole polarizability cannot be obtained from finite-field calculations, so it 
was obtained from the sum-over-state expression with the transition moments and excitation 
energies computed with the MRCI method. All calculations of the long-range coefficients 
employed both dalton 24] and molpro [25] suites of codes. 

The transitions from the ground X X E state to the X E and 1 II states are electric dipole 
allowed. The transition dipole moments for the electric, /ij, transitions were computed from 



the following expression 



27j: 

A*i=(*^ B W«SB> ) (10) 



where denotes the ith component of the position vector , while an d ^ab are the wave 
functions for the ground and excited states, respectively. Note that in Eq. ffTUl) i = x or y 
corresponds to transitions to 1 I1 states, while i = z corresponds to transitions to states. 
In the present calculations the electric transition dipole moments were computed as the first 



residue of the LRCCSD linear response function with two electric, r, operators [2JJ. In these 
calculations we have used the dalton program j^J]. We have evaluated the dependence of 
the transition dipole moments with the internuclear distance for the same set of distances 
as the excited state potential energy curves. 

As will be shown in the next section the electronic states of the low lying excited states 
of the (BaRb) + molecular ion show strong nonadiabatic couplings. Therefore, in this work 
we have computed the most important radial 

R(R) = ((n) 2S+1 \A\\^\(nf s+i \A\), (11) 



and angular 

A(R) = ((n) 2S+1 \A\\L+\(n') 2S+1 \A'\), (12) 

coupling matrix elements. In the above equations d/dR and L + denote differentiation 
with respect to the internuclear distance and the electronic angular momentum operator, 
respectively. Note that the radial operator couples states of the same multiplicity and 
symmetry, while the electronic angular momentum operator couples states with A differing 
by one. In the present calculations the angular coupling between the singlet states was 
computed as the first residue of the LRCCSD linear response function with two angular 
momentum operators L 21[. In these calculations we have used the dalton program 24]. 
All other nonadiabatic couplings were obtained with the MRCI method and the molpro 
code [25[. We have evaluated the dependence of the nonadiabatic coupling matrix elements 
with the internuclear distance for the same set of distances as the excited state potential 
energy curves. 

Barium and rubidium are heavy systems, so the electronic states of the (BaRb) + molec- 
ular ion are strongly mixed by the spin-orbit (SO) interactions. Therefore, in any analysis 
of the collisional cross sections between Ba + and Rb the SO coupling and its dependence on 
the internuclear distance R must be taken into account. We have evaluated the spin-orbit 
coupling matrix elements for the lowest dimer states that couple to the + and 1 states of 
(BaRb)" 1 ", with the spin-orbit coupling operator if so defined within the Breit-Pauli approxi- 
mation 36j . The spin-orbit coupling matrix elements have been computed within the MRCI 
framework with the molpro code [25j. The full spin-orbit Hamiltonian has been used in the 
calculations, i.e. both the one- and two-electron spin-orbit integrals were included. Having 
the spin-orbit coupling matrix elements at hand, we can build up the matrices that will 
generate the potential energies of the spin-orbit states that couple to + and 1 symmetry. 
The matrix for the 1 states writes: 



/ y(D ln ((l^niHsoKi)^) ((i)!n[if so |(i) 3 A) ((i) 1 n|H so |(2) 3 s) <(i) 1 n| J ff so |(i) 3 n) \ 

((ij^ihsoKi) 1 ^ y(D 3s o o <(i) 3 E|H so |(i) 3 n> 

{(^SAIHsoKl) 1 ^ y« 3A - <(l) 3 A|ff so |(l) 3 A) <(l) 3 A|H so |(l) 3 n} 

{(2) 3 E|H so |(i) 1 n> o o y(2) 3s <(2) 3 £|H so |(i) 3 n> 

V ((l^ni^soKi) 1 !!) <(i) 3 n|H so |(i) 3 £) <(i) 3 n||H so |(i) 3 A> <(i) 3 n|H so |(2) 3 s> yd) 3n 



(13) 
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(14) 



while the matrix for the + states is given by: 

( y(i) 3 n_((i)3n|tf so |(i)3n) ((i^niijsoix 1 !]) ((i) 3 n|i7 so |(2) 1 E) ((l^niffsoKi)^) ^ 

(XiEl^soKl) 3 ^ V x ^ 

((2) 1 E|ff S o|(l) 3 n) 

\ <(3) 1 S|i^so|(l) 3 n) ] 

Diagonalization of these matrices gives the spin-orbit coupled potential energy curves for 
the 1 and + states, respectively. Note that all potentials in the matrices (|T3|) and ( !T4|) 
were taken from CCSD(T), RCCSD(T), and LRCC3 calculations. Only the diagonal and 
nondiagonal spin-orbit coupling matrix elements were obtained with the MRCI method. 
Once the eigenvectors of these matrices are available, one can easily get the electric dipole 
transition moments and the nonadiabatic coupling matrix elements between the relativistic 
states. 

Finally, the long-range coefficients corresponding to the relativistic potentials were ob- 
tained by diagonalizing the matrices of Eq. (|T3|) and with the potentials expanded 
according to Eqs. (E])-(j3]) and the spin-orbit coupling matrix elements fixed at the atomic 



values. Note that unlike in the case of resonant interactions between like atoms [43], the SO 
coupling does not change the leading power in the asymptotic expansion of the interaction 
energy, but only changes the numerical values of the coefficients. One should also note that 
the atomic SO coupling does not change in our model the long-range coefficients for the 
(l) 1 ^ and (2) X S states due to the different dissociation limits: Ba( 1 / 3 D)+Rb + ( 1 S) versus 
Ba+( 2 S)+Rb( 2 S). 

In order to mimic the scalar relativistic effects some electrons were described by pseu- 
dopotentials. For Ba and Rb we took the ECP46MDF and ECP28MDF [45j pseudopo- 
tentials, respectively, from the Stuttgart library. For both barium and rubidium we used 



the spdfg quality basis sets suggested in Refs. |44], |45| . The full basis of the dimer was used 
in the supermolecule calculations and the Boys and Bernardi scheme was used to correct for 
the basis-set superposition error j^ . 

It should be stressed at this point that the ab initio results reported in the present paper 
obtained by the ab initio methods described above will allow to perform dynamical calcula- 
tions of the cross sections in the nonadiabatic, multichannel regime, both in the Hund's case 
(a) (nonrelativistic states, SO and nonadiabatic couplings, and transition moments) and (c) 
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(relativistic potentials, nonadiabatic couplings, and transition moments), cf. Ref. [46 ] . 

Finally, to conclude this section we would like to emphasize that almost all ab initio re- 
sults were obtained with the most advanced size-consistent methods of quantum chemistry: 
CCSD(T), RCCSD(T), and LRCC3. Only the SO coupling matrix elements and nonadia- 
batic matrix elements were obtained with the MRCI method which is not size consistent. 
Fortunately enough, all of the couplings are important in the region of the curve crossings 
or avoided crossings and vanish at large distances, so the effect of the size-inconsistency of 
MRCI on our results should not be drammatic. 



B. Dynamical calculations 

In collisions of the Ba + ( 2 S) ion with Rb( 2 S) atom we have basically three types of pro- 
cesses: elastic scattering in the singlet and triplet potentials, spin-flip (spin-exchange) pro- 
cess, and the inelastic radiative charge transfer from the singlet and triplet manifolds of 
Ba + Rb to the ground state of BaRb + . In the present paper we restrict ourselves to single 
channel calculations. A more detailed multichannel treatment will be presented elsewhere. 

To compute the elastic cross sections we need to solve the radial Schrodinger equation 
for the relative motion of the Ba + ion and Rb atom at an energy E: 

subject to the following normalization conditions: 

POO 

/ *% lJ (R)* JSa j(R)dR = 6(E 1 -E2), (16) 
Jo 

where ^ej(R) is the scattering wave function, /i is the reduced mass of the (BaRb) + ion, 
and V(R) stands for the interaction potential of Ba + ( 2 S) with Rb( 2 S) in the singlet or triplet 
manifold. Note that the normalization condition of Eq. (IT6|) is equivalent to the following 
large- R behavior of the wave function ^ej(R), 

2fi \ 1/2 . /, _ Jtt 



where Sj(E) denote the phase shift corresponding to the J partial wave, and the wave vector 

h 2 k 2 

2;u 



k is given by the standard expression, E = Equation (TToT) subject to the normalization 
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condition (TIB"]) allows us to compute the cross sections for the elastic and spin-flip collisions 
from the standard expressions: 

A 00 A 00 

<^ = -U E( 2J + !) sin2 <^)> = £ E( 2 J + !) sin2 <^)> ( 18 ) 

j=o j=o 

. 00 

4-7T 

^ = p E( 2J + !) sin2 W(^) - #(*0) > ( 19 ) 
j=o 

where the superscripts "s" and "t" on a e \ and <5j pertain to the singlet and triplet potentials, 
respectively. Note that an exact description of the spin-flip process would require at least 



two coupled channels, so the expression ( fl9|) is only approximate [37]] . It is derived under 
assumption that the hyperfme splittings are much smaller than the collision energy. However, 
it was shown to work relatively well, even at low energies 



38 



Theoretical description of the charge transfer process between the atom and the ion is 
somewhat more ellaborate. To the first-order of perturbation theory the radiative charge 



transfer can be described by the following Fermi golden type expression 40]: 



Am 2 ft 

*«(E) = -^T A ( E )> (2°) 
where the Einstein coefficient A(E) is given by: 

4 3 00 roc 

A(E) =^WE E (2 J/ + 1 )[/ £ 3 Hj>\(*E'jMR)\*E»J>>)\ 2 de 

J'=0J"=J'±1 Jo 

+ J2 H A E v''J" - ET\^eoMR)\^j-)\% (21) 



where the primed and double primed quantities pertain to the excited and ground state 
potentials, respectively, fi(R) is the transition moment from the ground to the excited elec- 
tronic state, a = 1/137.035999679(94) is the fine structure constant, e is the electron charge, 
and the Hohn-London factor Hj> is equal to ( J' + 1)/ (2 J' + 1) for the P branch (J' = J" — 1) , 
and to J'/(2J' + 1) for the R branch (J' = J" + 1). Here, e stands for the energy difference 

e = E" -E' + Ajp, (22) 

where Ajp is the difference of the ionization potentials. The scattering wave functions 
appearing in the expression ( l2Tj) are solutions of Eq. (fl5|) . while the bound-state wave 
functions fulfill the following Schrodinger equation: 

d 2 2a rr/ . J"(J" + 1) 2nE v nj»\ 



dR 2 h 2 v ' R 2 h 2 
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subject to the following normalization conditions: 

POO 

/ ^ VlJ (R)V V2 j(R)dR = 5 Vl 
Jo 



t'2 ■ 



(24) 



where V(R) stands for the ground state potential of the (BaRb) + molecular ion. 



A significantly simpler approach proposed in Ref. [41] approximates the sum over all 
continuum and bound states in Eq. (121 p with a simple average of a space-varying Einstein 
coefficient over the initial scattering wave function ^ e'j 1 '- 



A(E) = ^T(2J' + 1)(*e>j>\A(R)\*e>j>), 



J'=0 



where 



A(R) 



a' 



3he 6 



(6V) 3 (R)» 2 (R), 



(25) 



(26) 



and SV(R) is the difference between the excited and ground state potentials. 



III. NUMERICAL RESULTS AND DISCUSSION 

A. Nonrelativistic potential energy curves and spectroscopic characteristics of the 
ground and excited states 

Calculations were done for the ground state and first eight (four singlet and four triplet) 
excited states of (BaRb) + . Two states dissociate into Ba + ( 2 S)+Rb( 2 S), three states into 
Ba( 3 D) + Rb + ( 1 S), and three states into Ba( x D) + Rb + ( 1 S). The potential energies were 
calculated for twenty interatomic distances R ranging from 4 to 50 bohr for each potential 
curve. The potential curves are plotted in Fig. (TJ while the spectroscopic characteristics 
of these states are reported in Table fl] The ground state is absent on figures due to its 
regular behaviour and single minimum, and in order to increase the visibility of the other 
states. The separated atoms energy for each state was set equal to the experimental value. 
Numerical values of the potentials are available from the authors on request. 

Before going on with the discussion of the potentials let us note that the atomic excitation 
energies obtained from the LRCC3 calculations are very accurate. Our predicted position of 
the nonrelativistic 3 D state of barium is 9422 cm -1 , to be compared with the experimental 



value of 9357 cm 



-i 



421 ] deduced from the positions of the states in the D multiplet and the 



Lande rule. For the 1 D state of Ba we obtain 11907 cm in a relatively good agreement 
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with the experimental value of 11395 cm" 1 42]. It is worth noting that the present results 
for the atomic excitation energies are as accurate as the results of fully relativistic atomic 
calculations of Kozlov and Porsev [52], and more accurate than the data obtained from fully 



relativistic Dirac-Coulomb calculations 



201 ] . To further assess the quality of the methods, 



basis-sets and pseudopotentials employed in the present paper we have computed the static 
polarizabilities of the ground state of Ba atom, of the ground state of the Ba + ion, and 
the scalar and tensor components of the polarizability of Ba( 1 D). The present polarizabiity 
of the ground state of the barium atom is 272.5 a.u. The experimental value is 268 ± 22 



52] is 264 a.u.. Also the static 



, while the best theoretical result of Kozlov and Porsev 
polarizability of the Ba + ion, 132 a.u., is in a fairly goo d ag reement with the result of fully 
relativistic calculations of Ref. |49J, 124.15 a.u., of Ref. 50j, 124.26 a.u., and with the most 
recent experimental result, 123.88(5) a.u. [51]. For the X D state we get the scalar and tensor 



polarizabilities of 289 and 73 a.u., in a fair agreement with the results of Ref. [52], 266 and 81 
a.u., respectively. The methods employed in the present paper do not allow for a consistent 
calculation of the dissociation limit Aip of the (2) X S and (1) 3 S states corresponding to 
Ba + ( 2 S)+Rb( 2 S) separated atoms, since the ground state calculations employed CCSD(T) 
and the calculations on the singlet and triplet excited states the LRCC3 and RCCSD(T) 
methods, respectively. However, we can estimate Aip from the energy of the (2) 1 S and (1) 3 S 
at R = 50 bohr. Both the LRCC3 calculation on the singlet state and RCCSD(T) calculation 
on the triplet state give Arp = 8097 cm -1 in a good agreement with the experimental value 
of 8344 cm -1 42|. The fully relativistic result of Ref. 20], 8065 cm -1 , is very close to our 
value. 

The ground X X E state of the (BaRb) + molecular ion is a strongly bound state with the 
binding energy of 5136 cm -1 . The minimum on the potential energy curve for this state 
appears at a relatively large distance R e = 8.67 bohr. The origins of the binding can 
be explained by using the symmetry-adapted perturbation theory of intermolecular forces 
(SAPT) 28|, |29j]. As could be guessed, the interaction energy at the minimum results from a 
subtle balance of the induction attraction and exchange-repulsion. The induction energy is 
huge, -14499 cm -1 , but is strongly quenched by the exchange-repulsion term, 10068 cm -1 . 
The electrostatic contribution due to the charge overlap of the unperturbed electron clouds 
of Ba and Rb + , -554 cm -1 , and dispersion term, -598 cm -1 is of minor importance, again 
in agreement with our intuition. 
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An inspection of Fig. [T] shows that the potential energy curves for the excited states 
of the (BaRb) + molecular ion are smooth with well defined minima. The potential energy 
curves of the (2) and (S) 1 ^ show an avoided crossing. The potential energy curves of the 
(2) X S and (2) 3 E states exhibit a double minimum structure. The double minimum on 
the potential energy curve of the (2) X S state is due to the interaction with (3) 1 !] state. 
Other double minimum structure can be explained from the long-range theory and will be 
discussed below. Some potential energy curves show maxima. These are due to the first- 
order electrostatic interactions in the long range, and will also be discussed in more details 
in sec. IIIIBI Except for the shallow double minima structure of the (2) X S and (2) 3 S states, 
and a shallow (S) 1 ^ state, all other excited states of the (BaRb) + molecular ion are strongly 
bound with binding energies D e ranging from approximately 4380 cm -1 for the (1) 3 A state 
up to 6301 cm -1 for the (l) 3 n state. The (2) X S state, important from the experimental 
point of view, has two minima at R e = 9.02 and 15.19 bohr. The depths of these minima 
are 911 and 576 cm -1 , while the barrier separating them is of 30 cm -1 suggesting that the 
tunneling between the two wells will be very fast. The (2) 3 S state also shows two minima 
at R e 9.82 and 16.78 bohr of 1874 and 697 cm -1 , respectively, separated by a barrier of 681 
cm -1 . 

Before comparing our results with the ab initio data reported in Refs. [igL I20I ] let us stress 
the importance of the triple excitations in the wave functions for some states. This point is 
illustrated in Fig. El where LRCCSD versus LRCC3 and CCSD versus CCSD(T) potential 
energy curves are plotted for selected states. An inspection of these figures shows that the 
contribution of the triple excitations is important for the excited (2) 1 E and (2) 3 E states, and 
relatively unimportant for the (1) 3 A and for the (l) 3 n states. Not shown on these figures are 
the (1) 3 E, (l) x n, and (1) X A states. For these states we find that the potential for the (1) 3 S 
state is relatively unaffected by the triple excitations, while the two other potential energy 
curves show large differences depending whether the T3 cluster operator is included or not 
in the wave functions. These results strongly suggest that the CCSD or LRCCSD method 
works for those states that can be described by a single reference determinant. For open- 
shell states, i.e. all singlet states and the (2) 3 S state, the effect of the triple excitations 
is large and changes both the well depths and the barriers. In principle, the differences 
could be due to CCSD(T) vs. LRCC3 methods used in the calculations. However, as we 
already stated in sec. [Ill we have checked for a few points that the RCCSD(T) and LRCC3 



15 



results for the (1) 3 £, (l) 3 n, and (1) 3 A are very close. To get a better understanding of the 
importance of the triple excitations we analysed the energy gaps between the lowest electron 
configuration for a given state and the lowest triply excited configuration. It turned out that 
for the effectively high-spin states this energy difference was high, while for the manifestly 
states the opposite was true. Using simple perturbation theory arguments one can deduce 
that the correlation due to the triple excitations will be important for states with a small 
energy gap between the lowest and triply excited configurations, and much less important 
for states with large energy difference. 

Let us compare our results with other available ab initio data 19j. The spectroscopic 
isted in Table fl] and compared with the results of Ref. Unfortunately, 



constants are 
Knecht et al. 



19| did not report the binding energies D e of the molecular states, but only 



the values of the electronic term values T e taking the minimum of the ground state as zero 
of energy. An inspection of Table |T] shows that the agreement with the data of Knecht et al. 
191 ] is relatively good, given the fact that their results were obtained using the internally 
contracted multireference configuration singles and doubles method based on a CASSCF 
reference function. For most of the states the computed electronic terms agree within a few 
hundred cm -1 . For the high-spin states the differences in the positions of the minima are 
0.1 bohr at worst, while the electronic terms differ by 300 to 600 cm -1 . For the open-shell 
states, where the triple excitations in the wave function are important, the differences in 
the positions and well depths are more important. The most striking difference between 
the present results and the data of Ref. 19j is the (l) 1 !! state. Here, the difference in the 
position of the minimum is 0.6 bohr, and the difference in T e is as much as 3300 cm -1 . 
The double minimum structure of the (2) X E and (2) 3 S states was not reproduced by MRCI 
calculations, but the barriers and the avoided crossing between the (2) X S and (3) X E states 
are reproduced. In general, the agreement is only qualitative. To end this paragraph we 
would like to say that the results of fully relativistic Dirac-Coulomb calculations published 
by the same Authors in Ref. 20j are in a much better agreement with the present data, cf. 
sec. IfflCl 
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B. Nonadiabatic coupling matrix elements and electric transition dipole moments 
from the ground X X S state 

The most important nonadiabatic coupling matrix elements between the excited states 
of the (BaRb) + molecular ion are reported on Fig. [3J The regions where these couplings 
could possibly be important are indicated on Fig. [TJ The left pannel on this figure shows 
the radial coupling between the (2)*E and (S) 1 ^ states. The potential energy curves for 
these states reveal an avoided crossing shown on Fig. [1] at R ~ 12 bohr. An inspection of 
Fig. El shows that the maximum of the radial coupling between these two states corresponds 
to this distance. In general, the radial coupling as a function of the distance R is small, 
and rather localized around the point of the avoided crossing. The angular coupling matrix 
elements reported on the right pannel of Fig. H] show more variations with R. The angular 
coupling between the (l) 1 !! and (2) 1 S states has a broad maximum around R ~ 11 bohr, 
and this distance roughly corresponds to the crossing of the (l) 1 !! and (2) 1 E potential energy 
curves. The angular coupling between the (l) 3 n and (1) 3 S states shows a broad minimum 
at R ~ 7 bohr, and again this distance roughly corresponds to the crossing of the (l) 3 n and 
(1) 3 E potential energy curves. The last angular coupling that may influence the dynamics 
of the (BaRb) + molecular ion corresponds to the (2) 3 S and (l) 3 n states. Here, the R 
dependence of the angular coupling is quite different, but the largest variations correspond 
again to the point, where the two curves cross, R ~ 12 bohr. Note parenthetically that 
radial coupling tends to zero as R~ 7 , the angular couplings between the (2) 3 S and (l) 3 n 
states tends to a constant value, ( 3 D(Ml = 1)|L + | 3 D(M^ = 0)) = v3, while the coupling 
between the (l) 1 !! and (2) 1 S decays exponentially at large R. This unusual exponential 
decay is due to the different dissociation limits of the ground X X E and excited (2) X S states: 
Ba( 1 S)+Rb+( 1 S) versus Ba+( 2 S)+Rb( 2 S). It is gratifying to note that the MRCI method 
used in the calculations of the radial coupling and angular couplings between the triplet 
states quite precisely located the regions of the avoided crossing and curve crossing despite 
the fact the latter were determined from the RCCSD(T) and LRCC3 calculations. This 
suggests that the computed nonadiabatic coupling matrix elements are reliable, at least 
around of the crossings. 

The electric dipole transition moments between the ground state and the three excited 
electric dipole-allowed states, two and one 1 I1, are plotted in Fig. H]as functions of the 
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interatomic distance R. The calculated electric transition moments show a strong depen- 
dence on the internuclear distance R. For the transition moments to the excited states of 
the S symmetry the curves show broad maxima around the positions of the depths on the 
potential energy curves. The transition moment to the II state is very small, suggesting 
that this state will be of minor importance in the dynamics of the (BaRb) + ion. At large 
interatomic distances the transition moments to the (S) 1 ^ and (l) 1 !! states tend to zero 
as /i4i? -4 , while the transition moment to the (2) X S state decays exponentially with the 
internuclear distance R. This R~ 4 dependence can be derived fron the multipole expansion 
of the wave functions of the (3) 1 E, (l) 1 !!, and X X S states. The expressions for the leading 
long-range coefficient fi^ of the electric transition dipole moments read: 

/ i 4 = A| A |a Rb+ ( 1 S||g 2 || 1 D), (27) 

where A = 3/2^ and A 1 = A /\^3. 



C. Spin-orbit coupling and relativistic potential energy curves, nonadiabatc cou- 
pling, and electric transition dipole moments 

A large number of spin-orbit interactions couple the dimer states of the (BaRb) + molec- 
ular ion. In Fig. [5] we report the R dependence of the spin-orbit coupling matrix elements 
that couple to the 1 and + states. These states are most interesting for the collisional 
dynamics of (BaRb) + . Similar results can easily be obtained for the 0~, 2, and 3 states. 
An inspection of Fig. [5] shows that the SO coupling matrix elements have relatively large 
variations at small internuclear distances and tend to zero or to the atomic values at large 
R. All diagonal matrix elements of the spin-orbit Hamiltonian are lower than the atomic 
1 3 D spin-orbit constant of barium (« 150 cm -1 ), while the largest nondiagonal couplings are 
observed for the pairs of states associated with the crossing of the corresponding potential 
energy curves. The accuracy of the atomic SO couplings can be judged by comparison of the 
computed and observed positions of the energy levels in the 3 D multiplet. The calculated 
energies of the 3 Di, 3 D2, and 3 D3 states are 9035, 9254, and 9680 cm -1 , and are in a very 
good agreement with the experimental values of 9034, 9216, and 9597 cm -1 {42]. It is worth 
noting that some of the couplings vanish at large distances due to the different dissocia- 
tion limits: Ba+Rb + versus Ba + +Rb. This means that the neglect of the R dependence 
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of the spin-orbit matrix elements would lead to wrong relativistic potentials since some of 
important couplings would be neglected. As an example we can cite the coupling of the 
(1) 3 E and (l) 3 n states. The asymptotic value is zero due to the different dissociation limits 
of these states: Ba( 3 D)+Rb + ( 1 S) versus Ba + ( 2 S)+Rb( 2 S). Therefore, when approximating 
in Eq. (I13p all SO couplings by the corresponding atomic values one would obtain a com- 
pletely wrong matrix of the spin-orbit Hamiltonian for the 1 states since in the atomic limit 
((l) 3 S|iJso|(l) 3 n)=0. In particular, in the atomic approximation the (1) 3 S state would 
remain unchanged, and as it will be shown below it changes quite a lot at small internuclear 
distances. The same is true for the (2) 1 S state. 

The diagonalization of the spin-orbit Hamiltonian matrices, Eqs. ( I13p and (j!4j) . gives the 
potential energy curves of the states that couple to 1 or + symmetry. The corresponding 
curves for the + and 1 symmetries are reported on Fig. [61 An inspection of Figs. [1] and [6] 
shows that the crossings of the diabatic (nonrelativistic) states are transformed into avoided 
crossings on the spin-orbit coupled relativistic curves. Obviously, the inclusion of the spin- 
orbit interaction results in different dissociation pathways. Due to the presence of many 
closely located molecular states in the 3 D — X D energy range that couple to the + and 1 
symmetries, the effect of the spin-orbit coupling is very pronounced. Indeed, comparison of 
Fig. [1] and Fig. [6] shows that the behavior of some relativistic curves is drastically modified 
compared to the nonrelativistic case. 

In Table [III we report the spectroscopic constants of the relativistic states and compare 
them with the available ah initio data of Ref . [2o| • An inspection of this Table shows that 
the agreement between the present results and the data of Ref. 20J] is excellent for the 
dissociations Ba( 1 S )+Rb + ( 1 S ) and Ba + (' 2 S 1 /^.)+Rb(' 2 S 1 / 2 ) ; Indeed, our result for the well 



depth of the ground state overestimates the data of Ref. [20j by only 1.6%. For the first 
excited states of + and 1 symmetry, our results underestimate the values of Knecht et al. 
20[ by 3.1% and 3.7%, respectively. For all states mentioned above, the positions of the 
minima in the two calculations agree within 0.1 bohr or better. It is gratifying to observe 
such an excellent agreement between two different sets of ab initio calculations performed 
with different methods, CCSD(T) and RCCSD(T) in the present work vs. MRCI in Ref. 
[20! ] . Such a good agreement was expected from the analysis of the nonrelativistic results, 
since the triple excitations are relatively unimportant for these states. It is also worth noting 
that the pseudopotentials and basis sets used in our calculations 44, 4^ do a very good job, 
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as compared to the fully relativistic Dirac-Coulomb calculations 20] . The comparison for 
higher excited states is less favorable. For the (3)1 state the agreement of D e is within 4%, 
and R e is shifted by 0.1 bohr. This good agreement is again not fortuitous, since the (3)1 
relativistic state is dominated by the nonrelativistic (1) 3 A component, and the latter is a 
high-spin state not very sensitive to triple excitations in the wave function. For other states 
the differences in the well depths are of the order of 8% to 13%, and mostly reflect the lack 



of triple excitations in the calculations of Ref. [20] . We would like to stress, however, that 
the overall pictures of the relativistic states in the present paper and in Ref. [201 agree quite 
well. 

As in the nonrelativistic case, in the relativistic picture states of the same symmetry do 
not cross, while states of different symmetries can cross. Given the complicated pattern of 
the molecular states, cf. Fig. [61 the knowledge of the nonadiabatic couplings is essential 
for the multichannel dynamics. The nonadiabatic couplings between the relativistic states 
as functions of the internuclear distance R are presented on the left pannel of Fig. [7j An 
inspection of this figure shows that some couplings are rather localized in space with sharp 
maxima or minima, and some other show broad structure. All this can be rationalized 
by looking at the predominant singlet or triplet character of the states involved. Since all 
these structures can be explained in such a way, we take the coupling ((5)1|L + |(3)0 + ) as 
an example. The (5)1 and (3)0 + show crossing around R ~ 12 bohr, and in this region the 
nonadiabatic coupling has a broad maximum. At these distances both states are primarily 
singlets with only a small admixture of some triplet states. At R ~ 15 bohr the (3)0 + state 
shows an avoided crossing with the (2)0 + state. In the nonrelativistic picture the (2)0 + state 
is mostly dominated by the (l) 3 n state, while the (3)0 + state by the (2) X E state. Thus, at 
distances larger than the avoided crossing the matrix element of L + between the (l) 3 n and 
(l) x n states is zero. 

In the relativistic case many transitions that were forbidden at the nonrelativistic level 
due to the different multiplicities of the states involved become allowed due to admixtures of 
singlets to triplet and vice versa, cf. the right pannel of Fig. UJ These additional transition 
moments are very small, showing that the relativistic states obtained by admixture of the sin- 
glet states to triplets are almost pure triplet states. The transition moment ((l)0 + |z|(4)0 + ) 
resembles very much ((l) 1 S|z|(3) 1 S), with little differences only at small internuclear dis- 
tances. More interesting are the transition moments ((l)0 + |z|(2)0 + ) and ((1)0 + |,2|(3)0 + ). If 
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one would take the sum of these two curves, the resulting curve would strongly ressemble 
the graph for ((1) 1 E|^|(2) 1 E) transition moment, cf. Fig. HI The explanation of this fact is 
very simple. At distances smaller than R pa 15 bohr the (3)0 + is predominantly a singlet 
state, while (4)0 + has mostly triplet character. The situation is opposite at R larger than 
15 bohr. This explains why we get two curves with a sharp decay to zero around 15 bohr. 

D. Long-range behavior of the nonrelativistic and relativistic potentials 

When describing cold collisions between atoms it is crucial to ensure the proper long- 
range asymptotics of the interaction potential. The long-range coefficients describing the 
asymptotics of the nonrelativistic potentials are reported in Table IIHI An inspection of this 
Table shows that dispersion contribution Cg lsp is modest, but not negligible, of the order of 
5% to 20%. The induction and dispersion coefficients are always positive, so they describe 
the attractive parts of the potentials. The electrostatic coefficients differ in sign depending 
on the state considered and are responsible for the appearance of barriers and long-range 
minima on the potential energy curves. This is illustrated on the left pannel of Fig. [SJ 
where it is shown how the long-range asymptotics nicely fits the ab initio points. It is worth 
noting that the quadrupole moments are extremely sensitive to the electronic correlation. 
For instance, the value of Cf st for the (3) 1 !] state is 1.16 at the LRCC3 level and 2.05 with 
the LRCCSD method. This means that the triple excitations diminish Cf st by as much 
as 44%. Surprisingly enough, such an effect is not observed for the hexadecapole moment. 
Here, the C 5 elst at the LRCC3 and LRCCSD differ by only 6%. 

The long-range coefficients describing the large R asymptotics of the relativistic states are 
presented in Table IIVI As discussed in IIII CI in the atomic limit there is no spin-orbit cou- 
pling between the ground X 1 S, (1) 3 S, and (2) X S states with other states, so the long-range 
coefficients of the (1)0 + , (1)1, and (2)0 + states remain unchanged in this approximation. 
An inspection of Table ITVl shows the SO coupling of the 3 D and X D states modestly affects 
the long-range coefficients in the Ba( 1 D2)+Rb + ( 1 So) dissociation limit. The SO coupling in 
the 3 D multiplet introduces larger changes. For instance, the C*5 lst coefficient for the (2)1 
state is zero. This is not fortuitous, but only reflects the fact that the hexadecapole moment 
of an atom in a 3 D! state is identically zero. Again, the signs of the electrostatic coefficients 
are responsible for the barriers and long-range minima, cf. the right pannel of Fig. [SJ 
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E. Elastic cross sections, spin-exchange, and radiative charge transfer 

Thus far we have discussed the results of the electronic structure calculations. Now, we 
turn to the problem of sympathetic cooling of cold barium ions, and present the results 
for the elastic, spin-flip, and charge transfer cross sections for both the nonrelativistic and 
relativistic potentials. In Fig. [9] we report the elastic cross sections in the singlet and triplet 
manifolds, the spin-flip cross section, and the charge transfer cross section from the (2) X S 
to ground X X E state, all calculated from the nonrelativistic potentials. An inspection of 
this figure shows that in the ultracold regime the elastic cross sections behave according to 
the Wigner's threshold law, and the value of the cross section extrapolated to zero energy 
agrees very well with the one determined from the s wave scattering length. At energies 
around 100 nK the energy dependence of the cross sections starts to exhibit some structures 
related to shape resonances appearing due to the contributions of higher partial waves in 
the expansions ( jl8j) and to glory interference effects. Note that in the range of energies 
up to 1 niK the curves representing <jI\(E) and o^(E) are hardly distinguishable, despite 
the fact that the potential energy curves for the (2) X E and (1) 3 S states are quite different. 
This behavior is purely fortuitous, and is due to the fact that these two states have the 
same asymptotics and very close scattering lengths, a s and a t , equal to —3.53 x 10 5 and 
—4.26 x 10 5 A, respectively. The superelastic spin-flip cross section shows a qualitatively 
similar behavior. Overall, all the elastic cross sections are very large, from around 10 10 A 2 
at ultralow temperatures to ^10 6 A 2 in the mK region. 

In the ultacold regime the charge transfer cross section, which is an inelastic state chang- 
ing cross section, decays as i? -1 / 2 , in accordance with the Wigner's threshold law. At 
nanoKelvin temperatures this cross section, of the order of 10 4 A 2 is five orders of magni- 
tude smaller than the elastic cross sections. When we go up to the milliKelvin temperatures 
this ratio is even slightly more favorable, about five orders of magnitude of difference. Thus 
we can conclude that at the nonrelativistic level and with the single channel description 
of the collisional dynamics, cooling of the barium ion by collisions with ultracold rubidium 
atoms should be very efficient. 

The results of dynamical calculations on the relativistic (2)0 + and (1)1 potentials are 
presented on Fig. [TUJ The elastic cross section from the (1)1 state is almost indentical to 
the cross section obtained with the nonrelativistic (1) 3 S potential. This is not surprising 
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since the spin-orbit coupling has a very small effect on the potential, slightly moving the 
repulsive wall, cf. Figs. [1] and [6j The relativistic scattering length is lower compared 
to the nonrelativistic (1) 3 £ value, but most of the features, resonance structure and glory 
interference effects, are almost the same. By contrast, the energy dependence of the elastic 
cross section in the (2)0 + potential is very different from that presented in Fig. H] for 
the nonrelativistic potential. We note that the Wigner's limit is very different, and the 
resonant structure, glory ondulation, and all details at higher energies changed drastically. 
For instance, the scattering lengths for the (2)0 + and (1)1 states are of —1.58 x 10 3 and 
—2.56 x 10 5 A, to be compared with the values of a s and a t , —4.26 x 10 5 and —3.53 x 10 5 A, 
quoted above. All these difference are not surprising, however, since the spin-orbit mixing 
has a profound effect in this case. In fact, up to R ~ 15 bohr, the point of the avoided crossing 
between the potential energy curves of the (2)0 + and (3)0 + states, the (2)0 + potential has 
predominantly the character of the (1) 3 IT state. Only after the avoided crossing and mostly 
in the long range it becomes an almost pure (2) X S state. Also the spin-flip cross section 
computed from the relativistic potentials is quite different from the nonrelativistic one. It 
shows a sharp resonant structure around in the /iK region. 

The stricking difference between the nonrelativistic and relativistic pictures is the fact 
that the charge transfer process is now allowed from both the (2)0 + and (1)1 states. The 
corresponding cross sections as functions of the energy are also shown on Fig. |9j Again, 
we observe that the Wigner's threshold law with the E' 1 ^ 2 decay of the cross section in the 
ultracold regime is preserved. An inspection of this figure shows that the charge transfer 
process from the (1)1 state to the ground state will be very slow, unimportant at temper- 
atures interesting from the experimental point of view. This is not a surprise, since the 
transition dipole moment from the (1)1 state to the ground state is very small, cf. Fig. [7J 
The charge transfer cross section from the (2)0 + state is significantly larger, but again is 
several orders of magnitude smaller than the elastic one over all the range of temperatures 
considered in our work. In particular, in the milliKelvin regime, the inelastic events are 
six orders of magnitude less probable than the elastic. Thus, the relativistic calculations 
confirm the conclusions from the nonrelativistic case that the sympathetic cooling of the 
barium ion by collisions with ultracold rubidium atoms will be very efficient. 

We would like to conclude this section by saying that the present single channel cal- 
culations strongly suggest that the sympathetic cooling of the barium atom by collisions 
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with ultracold rubidium atoms should be very efficient in the temperature range up to mil- 
liKelvin. The analysis presented here neglects the effect of the hyperfine splittings which are 
important at low collision energies. Thus, our single channel analysis performed in terms of 
singlet and triplet properties is only approximate at ultracold temperatures, and the actual 
collision rates will depend on the hyperfine spin states of the atom and ion. Moreover, other 
inelastic processes due to the presence of other channels could enter the game, and final 
conclusions can be presented only when full multichannel calculations are performed. Work 
in direction is in progress [53 ]. 



IV. SUMMARY AND CONCLUSIONS 



In this paper we have reported theoretical prospects for the symphatetic cooling of the 
barium ion by collisions with the ultracold buffer gas of rubidium atoms. Potential en- 
ergy curves for the ground X X S state of the (BaRb) + molecular ion corresponding to the 
Rb + ( 1 S)+Ba( 1 S) dissociation and for the excited states, (1) 3 S and (2) X E, corresponding to 
the Rb( 2 S)+Ba+( 2 S) dissociation were computed by means of size-consistent coupled cluster 
methods with single, double, and triple excitations in the wave function. The inclusion of 
the triple excitations in the wave function was shown to be essential for the open-shell (2) X S 
state. The (BaRb) 4 " molecular ion shows a lot of low-lying molecular states, and the cor- 
responding potential energy curves were computed as well. Using the molecular spin-orbit 
coupling matrix elements relativistic potential energy curves were obtained. The asymp- 
totics of the nonrelativistic and relativistic potentials was fixed with long-range coefficients 
calculated from C3 up to and including Cq. A good understanding of the dynamics in 
the (BaRb) + system requires the knowledge of the nonadiabatic (radial and angular) cou- 
pling matrix elements, which were calculated as well. Finally, the electric dipole transition 
moments from the ground state were computed. Based on the above ab initio electronic 
structure calculations, single channel dynamical calculations of the elastic, spin-exchange, 
and inelastic cross sections for the collisions of the Ba + ( 2 S) ion with the Rb( 2 S) in the energy 
range from to 1 mK were performed, using both the nonrelativistic and relativistic poten- 
tials. It was found that the elastic processes are a few orders of magnitude more favorable 
that the inelastic ones. Thus, we can conclude our paper by saying that the sympathetic 
cooling of the barium ion by the buffer gas of ultracold rubidium atoms should be very 
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efficient taking into account the two-body collisional properties of Rb and Ba + . 
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TABLE I: Spectroscopic characteristics of the nonrelativistic electronic states of the (BaRb) 
molecular ion. 
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TABLE II: Spectroscopic characteristics of the relativistic electronic states of the (BaRb) + molec- 
ular ion. 
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TABLE III: Long-range coefficients (in atomic units) for the nonrelativistic electronic states of the 
(BaRb)+ molecular ion. C 6 is the sum Q nd + Cf sp . 
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TABLE IV: Long-range coefficients (in atomic units) for the relativistic electronic states of the 
(BaRb) + molecular ion. 



State 
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FIG. 1: Nonrelativistic potential energy curves for the excited states of the (BaRb) + molecular 
ion. 
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FIG. 2: Comparison of the LRCC3 with CCSD, and RCCSD(T) with LRCCSD potential energy 
curves for selected states of the (BaRb) + molecular ion. 
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FIG. 3: Nonadiabatic coupling matrix elements of the nonrelativistic electronic states of the 
(BaRb) + molecular ion. The left and right pannels correspond to the radial and angular cou- 
plings, respectively. 
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FIG. 4: Electric transition dipole moments from the ground X 1 !] state to the 1 S and 1 n states of 
(BaRb)+. 
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FIG. 5: Matrix elements of the spin-orbit coupling for the electronic states of (BaRb) + . The left 
and right pannels correspond to + and 1 states, respectively. 
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FIG. 6: Relativistic potential energy curves for the excited states of the (BaRb) + molecular ion. 
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FIG. 7: Nonadiabatic coupling matrix elements of the relativistic electronic states and relativistic 
electric transition dipole moments of the (BaRb) + molecular ion. The left and right pannels 
correspond to the nonadiabatic couplings and transition moments, respectively. 




Ba('D)+Rb*('S) 



Ba('D)H-Rb*('S) 



Ba*( 2 S)+RbfS) 



I — 1 — I — 1 — I — ' — I — ' — I — 1 — I — 1 — I 

1 5 20 25 30 35 40 45 50 

R (bohr) 



12000 
11000 
10000 
9000 
8000 
7000 
6000 
5000 
4000 
3000 
2000 



*"'-*-^Ba( , D,)+Rb fSJ 



* * ★ ! B. 



il'DjjtRb'l's,) 
|( 3 D 2 )+Rb + ('S ) 
■('D.l+RbYS,) 
i*fs„J*Rb( ! S„J 



— i — 1 — i — 1 — i — ' — i — 1 — i — 1 — i — 1 — r~ 
15 20 25 30 35 40 45 50 

R (bohr) 



FIG. 8: Long range nonrelativistic (left pannel) and relativistic (right pannel) potential energy 
curves of the (BaRb) + molecular ion. Stars are based on long range coefficients like in Eq. [3j 
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FIG. 9: Elastic, spin-flip, and charge transfer cross sections for collisions of 138 Ba + ( 2 S) and 
87 Rb( 2 S) as functions of the collision energy from nonrelativistic potentials. 
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FIG. 10: Elastic, spin-flip, and charge transfer cross sections for collisions of 138 Ba + ( 2 S) and 
87 Rb( 2 S) as functions of the collision energy from relativistic potentials. 
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